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I. INTRODUCTION 

The problem of quantizing gravity has been debated frequently during the past decades. In doing so, it was 
established that the quantization of Einstein’s general relativity inevitably results in fundamental problems such as 
the perturbative nonrenormalizability [l], This represents a motivation to introduce other, more radical approaches 
to obtain a quantum theory of gravity, including the consideration of higher-order theories of gravity, string theory, 
and loop quantum gravity (2j. 

One possibility is to reject the perturbative quantization, as is done, for example, within the framework of loop- 
space nonperturbative quantum gravity d]. On the other hand, in order to quantize gravity one can try to adopt the 
nonperturbative quantization technique employed by Heisenberg when considering a nonlinear spinor field theory d ■ 
The central idea of this approach is that the description of the quantum system is achieved by using an infinite set of 
equations for all Green’s functions. Based on physically motivated arguments, one can introduce a cutoffprocedure 
to obtain a finite number of equations. Such a procedure is close to the one used in turbulence modeling [5[. 

Working within this approach, we here consider the situation where a quantum metric can be decomposed into a 
sum of an averaged (classical) metric g^ and a fluctuating (quantum) part 8g^ v . Also, we assume that the Green’s 
function of a product of the fluctuating part of the metric can be approximated in a certain way. Under these 
assumptions we calculate the action for the gravitational field and for matter (in the form of a scalar field) with an 
accuracy of 8g 2 . We show that the resulting action represents modified F(RY- type theories of gravity, which are now 
widely used in modeling the present accelerated expansion of the Universe [6(. 

For the quantum gravitating physical system considered here, with the decomposition of the metric into classical 
and quantum parts, the expectation value of the quantum part can be taken to be zero or nonzero, based on certain 
physical motivations. In the present paper we take it to be zero. This is in contrast to our previous work 7j, where the 
decomposition into two parts (classical and quantum) was also performed, but the expectation value of the quantum 
part of the metric was taken to be nonzero. Therefore, we here take the second variation of the metric into account, 
as discussed below. 


II. NONPERTURBATIVE QUANTIZATION TECHNIQUE 

According to Heisenberg’s nonperturbative quantization technique the classical fields appearing in the corresponding 
field equations are replaced by operators of these fields. For general relativity, one then has the operator Einstein 
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equations 


G p i/ - RfllS 


8 ttG - 

4 J /«o 


(1) 


where all geometric operators R pl/ , R p apv , T p are defined in the same way as in the classical case, and differ only 
in the replacement of the classical quantities by the corresponding operators (for details, see Ref. Q). 

There are no known mathematical tools for solving the operator equation ©■ The only possibility to work with such 
an operator equation is to average the equation © over all possible products of the metric operators g{x i) • • • g{x n ), 
and thus obtain an infinite set of equations for all Green’s functions: 


87tG 
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g{x i) • G 


fils' 
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g{x 1 )g{x 2 ) ■ G 


Ills' 


the product of g at different points (aq, • • • , x n ) ■ G 
the product of g at different points (aq, • • • , x n ) ■ 



8 t rG 
c 4 
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9{xi) ■ T, 


fiv 
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g(x 1 )g{x 2 ) ■ f p 
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( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


where | Q) is a quantum state (for details see Refs. 0,Q). The exact definitions of a non-perturbative vacuum state 
| 0) and a quantum state | Q) are given in Appendix [B] 

Evidently Eqs. ©I-© cannot be solved analytically. Different possibilities to solve them approximately were 
discussed in Refs. [2H9(. Note, that similar mathematical problems appear in turbulence modeling, where an infinite 
set of equations for all cumulants also arises [f|. 

Here we use the following strategy for approximate solving Eqs. ©-©: we decompose the metric operator g^ 
into classical, g pl/ , and quantum, Sg piy , parts and evaluate the expectation value of the Lagrangian with an accuracy 

. This strategy is similar to the one employed in connection with quantum torsion in Ref. [jjj]. 



III. ASSUMPTIONS OF THE QUANTUM AVERAGING 


In accordance with the quantization procedure, a metric in quantum gravity is an operator g pu . Here we consider 
a system for which the following decomposition is approximately valid 


~ SV + Sg + S 2 g 


( 6 ) 


Sg 


fils' 


Q) = 0; (Q 


S 2 9uis Q) 7 ^ 0; | Q) is some quantum state; Sg 


fils' 


where g^ is the classical part of the metric; \ Q 
and 5 2 g lll/ are the first- and second-order deviations of the operator g fll ,. The expression © is a quantum variant of 
the decomposition of the classical metric (see Ref. 0 . p. 129), where 1 Sg„ v is the first order term and S 2 g„ v is the 


second order term. In our quantum case it means that 


(sy 


« 2 9>- 


In Ref. t 7l we have considered a physical quantum system with the decomposition g^ u = g lw + Sg jlv . where g M! , is the 
classical part and Sg^ v is the quantum part of the metric with the nonzero vacuum expectation value (^Sg^^j. Thus 
the difference between the physical system of Ref. ,7) and the one of the present paper is that in the first case the final 
effective Lagrangian is calculated with accuracy ( Sg^ v j ^ 0, and in the second case the final averaged Lagrangian is 


calculated with accuracy / (Sg^ \ and (s 2 g /JiU 


For our approximate nonperturbative calculations, we insert the decomposition © into the Einstein-Hilbert action 
and evaluate it with an accuracy ^(Sg'j ^ ~ ^ S 2 g To do this, we have to make some assumptions on the 2-point 
Green’s function. Namely, we suppose that it can be decomposed as the product 


G 2 

; fils', pa (xi,x 2 ) (Q (^ 1 ) ■ Sg pcj {x 2 ) 
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( 7 ) 
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It is seen from this expression that G 2 should be symmetric under the permutations [i o v , p o a , and /x, v o p, a. 
In the subsequent calculations we will use the following assumptions: 

• The Green’s function G 2 (xi, x 2 ) can be approximately decomposed as the product of some tensors at the points 
X\ and x 2 : 


G 2 ;fii',p(T (xi,X 2 ) ~ Ppv(xi)P pa (x 2 )- ( 8 ) 

• Taking into account the symmetry properties noted above, one can see that there exist the following possibilities 
for choosing the tensor P pv : 

— P pi/ is proportional to the metric tensor: 

Ppv cx cjpi/- (9) 

— P^v is proportional to the Ricci tensor: 

(io) 

• The proportionality coefficient in the expressions © and ms should be some invariant. Consequently, it has 
to have the form F(R, R iLV R! LV , • • •). The coefficient F should be very small as g pv —> r where 77 ^ is the 
Minkowski metric. 

• The expectation value ^ 5 2 g lll ^ is given by some tensor of rank two, 

(&g^)=K^. ( 11 ) 

• Analogously to da and da, we assume that 

— K pv can be proportional to the metric tensor: 

Kfj,v cx g^v- ( 12 ) 

— K pv can be proportional to the Ricci tensor: 

(13) 

• P/J.V and K^,, can be a linear combination of the metric and Ricci tensors. For example, they could be the 
Einstein or Schouten tensors. 

• For each quantum state \Q) there exists only a single set of functions F and K pv . 

Thus, we assume that the quantum correlation between fluctuations of the metric at two points can be approximately 
described as 


G 2; „ (xi,x 2 ) « [P^F(R, R^RT, ■ ■ ■ )] ai • [P P „F(R, R^R^, ■ ■ • )] Ia . 
The simplest choice for F is 


(14) 


F = F(R ), 


(15) 


which corresponds to F(R)-gravities. 
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IV. EVALUATION OF THE AVERAGED ACTION 


We start from the classical Einstein-Hilbert Lagrangian 


C = -2^^ R ' 


( 16 ) 


where x = 87 rG/c 2 . We then expand C(g + Sg + 5 2 g) into a Taylor series and subsequently replace the classical 
quantities Sg,5 2 g by quantum ones Sg,S 2 g 


£(g + $g + S 2 g) w £{g) + 


SC 


Sg» v 


S 2 C 


-Sgt J,v 5g f ‘ 


S 2 C 


S 2 gP u . 


(17) 


Sg^ a 5g^SgP a 5 2 gP 1 ' 

Our next step is to average this Lagrangian over quantum fluctuations of the metric Sg F . Since in the present paper 
we assume that (^Sg lu ^ = 0 (see the Introduction and cf. Ref. where it was taken to be nonzero), we have 


6 2 £{g)) = (Sg^Sgpt 


S 2 C 

SgP v SgP° / ' S 2 gP’ v 

Taking into account the first variation of the Einstein-Hilbert Lagrangian 


5 2 gP 


SC(g) = -^^g(R, v - \g^R ) 


we can calculate the second variation as follows: 


S 2 C(g) = 


2 x 


(dyj yRf_u/ 2 d^Rj ~b a/ ~9^ yRfu/ ^Qfii/R 


5g^ + V=9G, v S 2 g^, 


(18) 


(19) 


( 20 ) 


where G^ = R^ — ^g^R is the Einstein tensor. Equation (f2U|) can be rewritten as 


^fiv 2 


( 


2 {Rfiv 2 Qh-uR \ Qat-flSg ^ $Rf±v 2 2 9^vSR 


Sg^ + G^S 2 g^y ( 21 ) 


Now we can calculate an expectation value of the Lagrangian (ED by replacing all classical quantities Sg p v and 8 2 gP 1 ' 

_ fll/ LLV _ _ 

in Eq. ([21]) by the quantum ones, Sg and 5 2 g (for details see Appendix [A]) . For simplicity, let us here consider 
the Ansatz ©, for which we obtain 


C(g + 5g + 5 2 g)j » - — ^~g[R - 2 RF(R, • • •) + 3 F{R, ■ ■ ■ )V^V A1 F(R, • • •) + G^K^) . (22) 

Thus, we see that we have derived a modified gravity theory. For the simplest choice F(R, ■ ■■) = F(R) and K = 0 
we have F(R )-gravity theory. 

Let us now perform similar calculations for the matter Lagrangian, using the decomposition given by Eqs. © and 
©• For simplicity, consider the scalar field (ft with the Lagrange density 


Cm — ~gCm — V* ~g 


V^V M 0-U(0) 


Its first variation is 


SCm = Ai? [V^V^ - g^L m ] SgT = 


where is the energy-momentum tensor. Then the second variation yields 


S 2 C m - | ^ 


9fia9i//3 + 2 9/J.iy9a/3 ) A m 9ct(3 ^ n4^ v 4* 


5g^5g ap + T^S 2 g^ 


(23) 


(24) 


(25) 


Using Eqs. ©-© and ED and replacing again all classical quantities Sg 1 * 1 ' and S 2 gP l/ by the quantum ones, we find 
the expectation value (c m + S 2 C m \ in the form 


(26) 


Cm + S 2 C m > = J~g 11 V'VV^ - [1 + 2F (A, • • •)] V(4i) + K ^ 


Thus, we see that a nonminimal coupling between the scalar field and gravity appears. Notice also, that in the case 
of ^ 0 one can obtain a gravitational theory in which the derivative of the scalar field (f>, appearing in the term 
T M! ,, is nonminimally coupled to curvature (for cosmological models with such a coupling, see, e.g., Refs. EHH). 
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V. CONCLUSION 

In the present paper we have considered the case of a quantum gravitating system when the metric can be de¬ 
composed into classical and quantum parts. For such a system, we have calculated the gravitational and matter 
Lagrangians with an accuracy up to the second variation of the metric. In doing so, we have decomposed the operator 
of the metric into a sum of its expectation value (c-number) and deviations (g-numbers) from this expectation value. 

Based on certain assumptions on the dispersion of quantum fluctuations of the metric, we have shown that: 

• Einstein gravity is modified in the spirit of E(I?)-gravity theories. 

• Matter is nonminimally coupled to gravity. 

In obtaining these results, we have assumed that the expectation value of the product {Sg, Jtv {xx) 8 gp ll/ {x 2 )^ at two 
points X\,X 2 can be decomposed into the product of two factors of some tensor at these two points, see Eqs. 0 and 
©• 

The proposed model, which takes into account the quantum fluctuations of the metric, provides us with the following 
scheme for the explanation of the present acceleration of the Universe: the quantum metric —> E(I?)-gravity —> the 
accelerated Universe. The model can explain qualitatively the smallness of the effective A-term, which comes from 
small quantum fluctuations of the metric. 

Note that the proposed procedure of quantizing the metric may be considered as being related to quantum gravity 
like the phenomenological Ginzburg-Landau model of superconductivity is related to the microscopical Bardeen- 
Cooper-Schrieffer theory: the proportionality coefficient F(R,---) can be calculated only from the true quantum 
gravity. 
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Appendix A: Variation of geometrical quantities 


Following Ref. 0 , we here give all formulae concerning the variation of R , R^, and R pl , a p- The variations of the 
inverse metric and the Christoffel symbols T^ p are 


Sg afi = -g ap g^Sg p(J = -h a ^ 

^ a.hfi'y ^ P^ycx. ^ 'yh'afi') “1“ ^7 

Sn p = \ (V v Sg£ + V„Sg£ - WSg v p ) , 


(Al) 

(A2) 

(A3) 


where g^ v and g are used to lower and raise indices, and for brevity, we have introduced = Sg Then the 
variation of the Riemann and Ricci tensors and the curvature scalar are 

SR a M5 = V 7 AT^ - \7s6T^, (A4) 

SRafi'yS 7 ^ hapR T h( 3 p R a ryfi (V /nh^ a V \7 ry Vq/^ 5/3 V 7 VV( 5 VqA 7 / 3 ) ; (A5) 

8 R a p = — —V P V P h a p + R a a h a p + Rp a h aa — 2 R ap p a h pcr + V Q V p hp P + V p h a p — V &Vph , (A 6 ) 

SR = h a PR a p + g af> 5R af} = h ap R a0 - \7 a \7 a h + V“V / 3 /i a ^, (A7) 

where h = g a ^h a p and the covariant derivative V, t is taken with respect to the metric g. 
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Let us now calculate the expectation values for the case ©: 


-77 (- o9^ R ) 9<xp(fig^6g a0 ) = 2 rf 2 {R,---), 




-(SR^Sg^) = -3 FV^V P F, 


—R^g^Sg^) = —2RF 2 (R, ■ ■ ■), 
~g^(6R6^) = -2RF 2 (R,---)+6FV^V p F, 


where we have used the fact 


() $9 pa-} = lim (V M ) Xi ^ 8g a p(xi)6g pa (x2)} • 


Summing (1A81) - (1A11I) yields 


^c) = c -V=~gRF{R,---). 


For the case we obtain the following expectation values: 


- o ( R i‘- ~ ) 9«f) (Sg^Sg^) = ~ ( R aP R af) - -R 


F 2 

~R 


\{dR^5g^) = 


-F- 
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V P V, 


FR, 


a/3 
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j-p 2 

+ 2 R°R afi R<*P — - 
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j ?2 dcc/3 

^R apfi aR p ° R af * Jp + 2F—V„V 


OL v P 
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R Q P 

F—\7 a VpF 

K 

1 „ / r— .— \ F 2 R a/ 3 R a/3 

"2 R \ 6 g^ V ) = 2 R ’ 

1 /—r—A F 2 R a pR a P . F__ F 


-2^\ SRS9 n = - 2 


R 


2 V “ V “ f - 2 


V q V /3 F^^ ) . 


Summing (1A14D - (1A17f) . wc have 




3 „n R n f<R a P 
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R - 
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1. 


-V Q Vp | F^4 | _ -V"V„ [ F 


F. 


-V^V ' 3 I F 


,R 


a/3 
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+ Ir-R ct/ 3 R“^ - Ir^R^R^ + f V'V P F + ^ 


( p R <±P^ 

-7 Vo-VflF 
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4 p 


}■ 


(AS) 

(A9) 

(A10) 

(All) 

(A12) 


(A13) 


(A14) 


(A15) 

(A16) 

(A17) 


(A18) 


Appendix B: Nonperturbative vacuum: discussion and definitions 

In perturbative quantum field theories a vacuum is defined by using an annihilation operator a as 

a| 0) = 0. (Bl) 

This definition explicitly uses the notion of quantum and consequently cannot be used for the definition of a nonper¬ 
turbative vacuum. Physically, the difference between perturbative and nonperturbative vacua is the following: the 
perturbative vacuum is a sea of virtual quanta that appear and annihilate everywhere and always; the nonperturbative 
vacuum is similar to a stormy sea with random waves (fluctuating fields) on it. 

We give the following definition of a nonperturbative vacuum for gravity (here we work with tetrads): 
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1. The expectation value of tetrad operators e° at any point x M should satisfy the following relation: 

(0 \e a/ j,(x p )e a l/ (x p )\ 0) = r)^, (B2) 

where rj^ u is the Minkowski metric. 

2. For some combinations of indices a, b , g,, ix the 2-point Green’s function of tetrad operators e° M is non-zero, 

G$ P> * {x p ,x°) = <0 |e° #t (*")e 6 I/ (^)| 0) ^ 0. (B3) 


3. The dispersion of the metric at any point x^ is non-zero, 


(e a/ x( x p )e a v (x p ) - T]^ 


0)^0. 


(B4) 


Thus, the definition of a nonperturbative vacuum for gravity is obtained by combining the set of equations for all 
Green’s functions (|B6j) - (|Bllj) with the constraints (1B2I) - (1B4I) . 

In order to obtain the set of equations for all Green’s functions, we use the operator Einstein equations 


If l V - mT^l,. 


The equations for all Green’s functions, which follow from the operator equation, are 


e ai ai K) 


e ai ai «)-e ^ 2 (x$)- 


e ai ai K)"-e a S n «) 


0y = x(^0 
0^ = 
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fLlS 
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0) = ( 0 


e ai ai ( x i )' (x%)-T liV 

e ai ai (^i)-‘-e a " n «)-T) 


(B5) 

(B6) 

(B7) 

(B8) 
(B9) 
(B10) 
(BH) 

with the constraints (IB2I) - (IB4I) . Equations (IB6MB11I) are partial differential equations for all Green’s functions. The 
solution of this set of equations with the constraints (1B2MB4I) gives us all Green’s functions describing a vacuum state 
in quantum gravity. The knowledge of all Green’s functions is identical to knowing the properties of the operators 
e a ll and the vacuum quantum state | 0). 

In the same way, we can define the quantum state | Q ): The quantum state | Q) and the properties of the operators 
e“ M are defined through all Green’s functions satisfying the set of equations 
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Qj = 

Q^j = x(q 

q\ = mIq 


T,. 


Q), 


Q), 


e ai ai K) • f M 


e ai ai ( x i) ■ e a i 2 O 2 ) 


Q), 


Q) = x\Q 


e ai ai K)"-eX { x n)' T ^ 


Q 


(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 


Let us emphasize that for the definition of the quantum state | Q) we do not use the constraints 

Finally, let us address the procedure of quantum averaging. In quantum mechanics, this procedure is carried out 
by the integration: 


L)= if* Lip dV. 


(BIB) 
























































For the non/perturbative quantization, this procedure is not obvious. To define a quantum average for an operator, or 
some product of operators, we first have to solve the set of equations (1B12I) - (1B17I) . Then, among the obtained Green’s 
functions, we find the required quantum average. 
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